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1.  INTRODUCTION 

The  stochastic  Petri  net  (spn)  model  is  well  suited  to  formal  representation  of 
concurrency,  synchronization,  and  communication  (cf.,  Marsan,  Conte,  and  Balbo  [11], 
Molloy  [14,15],  Natkin  [16],  Symons  [18]).  An  spn  is  specified  by  a  finite  set  of  places 
and  a  finite  number  of  transitions  along  with  a  normal  input  function ,  an  inhibitor  input 
function ,  and  an  output  function  (each  of  which  associates  a  set  of  places  with  a  transition). 
A  marking  of  an  spn  is  an  assignment  of  zero  or  more  tokens  to  the  places  in  the  net.  A 
transition  is  enabled  whenever  there  is  at  least  one  token  in  each  of  its  normal  input 
places  and  no  tokens  in  any  of  its  inhibitor  input  places;  otherwise,  it  is  disabled.  A 
transition  fires  by  removing  one  token  (per  place)  from  a  random  subset  of  its  normal 
input  places  and  depositing  one  token  (per  place)  in  a  random  subset  of  its  output  places. 

Heuristically,  an  spn  changes  marking  in  accordance  with  the  firing  of  a  transition 
enabled  in  the  current  marking.  Each  of  the  transitions  enabled  in  a  marking  competes  to 
change  the  marking  and  each  of  these  enabled  transitions  has  its  own  stochastic 
mechanism  for  determining  the  next  marking.  At  each  firing  of  a  transition  in  the  spn, 
new  transitions  may  become  enabled.  For  each  of  these  new  enabled  transitions,  a  clock 
indicating  the  time  until  the  transition  fires  is  set  according  to  an  independent  stochastic 
mechanism.  (There  is  no  restriction  to  exponentially  distributed  transition  firing  times.) 
If  an  enabled  transition  does  not  trigger  a  marking  change  but  is  enabled  in  the  next 
marking,  its  clock  continues  to  run;  if  such  a  transition  is  not  enabled  in  the  next 
marking,  its  clock  reading  is  abandoned. 
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The  spn  representation  provides  a  means  of  generating  sample  paths  for  the 
underlying  stochastic  process  of  a  discrete  event  simulation.  This  representation  is 
particularly  useful  in  connection  with  non*Markovian  systems.  The  "state  of  the  system 
at  time  /"  defines  the  underlying  stochastic  process  of  a  discrete  event  simulation.  When 
the  current  state  of  the  stochastic  process  is  s  and  a  previously  scheduled  event,  e*, 
occurs,  the  process  moves  to  a  new  state,  s.  The  marking  of  an  spn  corresponds  to  the 
state  of  the  process  and  the  firing  of  a  transition  corresponds  to  the  occurrence  of  an 
event.  The  graphical  representation  (bipartite  graph  of  places  and  transitions)  of  an  spn 
is  particularly  useful  in  that  it  incorporates  a  considerable  amount  of  information  about 
the  set  of  events  that  can  occur  when  the  process  is  in  state  s  and  the  sets  of  "new 
events"  and  "old  events"  when  event  e  triggers  a  transition  from  state  s  to  state  s. 

Although  steady  state  estimation  for  an  arbitrary  spn  is  a  very  difficult  problem, 
Haas  and  Shedler  [4]  have  provided  estimation  procedures  for  spn’s  that  are  regenerative 
processes.  To  establish  the  regenerative  property  for  an  spn,  it  is  necessary  to  show  the 
existence  of  an  infinite  sequence  of  random  time  points  at  which  the  process 
probabilistically  restarts.  It  is  often  clear  that  an  spn  probabilistically  restarts  when  a 
particular  transition  fires  leaving  the  system  with  a  fixed  marking.  For  specific  models, 
however,  it  is  nontrivial  to  determine  conditions  (distributional  assumptions)  under 
which  this  occurs  infinitely  often  with  probability  one.  Using  recurrence  theory  (Haas 
and  Shedler  [3])  for  generalized  semi-Markov  processes  (KOnig,  Matthes,  and 
Nawrotzki  [8,9],  Matthes  [12],  Whitt  [19]),  conditions  are  given  in  [4]  which  ensure  that 
an  spn  is  a  regenerative  process  in  continuous  time  with  finite  expected  time  between 


regeneration  points. 
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In  this  paper  we  focus  on  spn’s  with  special  structure  and  define  a  symmetric  spn. 
Informally,  an  spn  is  symmetric  if  there  are  mappings  of  places  onto  places,  markings 
onto  markings,  and  transitions  onto  transitions  which  preserve  the  sets  of  enabled 
transitions,  the  new  marking  probabilities,  the  sets  of  new  transitions,  and  the  clock 
setting  distributions.  Symmetric  spn’s  have  application  to  representation  of  ring 
networks  with  equally  spaced,  identical  ports;  cf.  Loucks,  Hamacher,  and  Preiss  [10]. 

Section  2  provides  the  formal  definition  of  an  spn  given  in  [4]  along  with  conditions 
which  ensure  that  an  spn  is  a  regenerative  process  and  that  the  expected  time  between 
regeneration  points  is  finite.  Using  a  geometric  trials  recurrence  criterion  (Iglehart  and 
Shedler  [6]),  Proposition  (2.19)  postulates  the  existence  of  a  transition,  e,  and  a 
marking,  jg,  such  that  transition  e  fires  and  the  new  marking  is  Jq  infinitely  often  with 
probability  one.  Conditions  on  the  old  clocks  ensure  that  the  process  probabilistically 
restarts  at  these  times.  This  result  is  the  basis  for  regenerative  simulation  of  spn’s. 

Section  3  considers  the  steady  state  estimation  problem  for  symmetric  spn’s.  Under 
the  assumptions  of  Proposition  (2.19),  regenerative  cycles  defined  by  the  times  at  which 
transition  e*  fires  and  the  new  marking  is  jg  can  be  decomposed  into  independent, 
nonidentically  distributed  blocks.  We  show  that  point  estimates  and  confidence  intervals 
for  characteristics  of  symmetric  functions  of  the  limiting  distribution  can  be  obtained  by 
simulating  the  symmetric  spn  in  blocks. 

In  Section  4  we  develop  estimation  procedures  for  passage  times  in  the  spn  setting. 
Formal  specification  of  a  sequence  { P'n:n  £  1}  of  passage  times  in  a  symmetric  spn 
{.XU):*  i  0}  with  marking  set,  S ,  and  transition  set,  £,  is  in  terms  of  four  subsets 
(^lti42,  2?j,  and  fl2)  The  sets  B1  and  #2  define  the  random  times  [7 jtj  i  1}  at 
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which  a  passage  time  terminates.  (The  sets  A1  and  A2  define  the  random  times  at  which  a 
passage  time  starts.)  Proposition  (4.1)  postulates  the  existence  of  e  eE  and  SqJqCS 
such  that  transition  e  fires  and  the  marking  changes  from  j0  to  Jq  infinitely  often  with 
probability  one  and  these  transition  firing  times  correspond  to  termination  of  a  passage 
time  with  no  other  passage  times  underway.  Conditions  on  the  "old  clocks"  ensure  that 
{(JK7i)tPj|+1):ii  2  0}  is  a  regenerative  process  in  discrete  time  and  that  the  expected  time 
between  regeneration  points  is  finite. 

Section  4  provides  two  estimation  procedures  for  passage  times  in  a  symmetric  spn. 
Each  of  these  procedures  rests  on  the  assumption  that  there  exist  e  eE,  sQe  Bv  and 
jq£2?2  satisfying  the  conditions  of  Proposition  (4.1).  The  regenerative  structure 
guarantees  that  Pn+P  as  n-+ao  and  the  goal  of  the  simulation  is  the  estimation  of 
r(J)  *  E{JXP)\,  where  / is  a  real-valued  measurable  function. 

Estimates  for  Kf)  can  be  based  on  measurement  of  passage  times  {P*:n  2:  1}  (a 
particular  random  subsequence  of  {P'n’-n  2  1})  and  simulation  of  the  underlying  spn  in 
regenerative  cycles.  Alternatively,  exploiting  properties  of  a  symmetric  spn,  estimates 
can  be  based  on  measurement  of  passage  times  2  1}  and  simulation  of  the 

underlying  spn  in  independent,  nonidentically  distributed  blocks.  This  estimation 
procedure  extracts  more  passage  time  information  from  a  simulation  of  fixed  length  and 
should  provide  estimates  for  r(f)  that  are  relatively  more  accurate.  In  Section  S  we 
verify  that  this  is  indeed  the  case  by  showing  that  the  resulting  confidence  intervals  are 


2.  REGENERATIVE  STOCHASTIC  PETRI  NETS 

Following  Haas  and  Shedler  [4],  formal  definition  of  an  spn  is  in  terms  of  a  general 
state  space  Markov  chain  (gssmc)  which  describes  the  process  at  successive  epochs  of 
transition  firing.  Let  D  =  {1,2,... ,7.}  be  the  index  set  for  a  finite  collection  of  places  and 
let  E  ■*  ^  a  finite  set  of  transitions.  Denote  by  S  the  finite  or  countable  set 

of  markings  and  for  se  S  write  s  =  (svs2,..., s^),  where  Sj  is  the  number  of  tokens  in  place 
j,  je  D.  Denote  the  index  set  of  the  normal  input  places  for  transition  eC  E  by  1(e)  £  D, 
the  index  set  of  the  inhibitor  input  places  by  L(e)  £  D,  and  the  index  set  of  the  output  places 
by  /(e)  £  D.  We  assume  that 

Ue)  fl  7(e)  =  0 

for  all  ee  E.  For  s  =  (sl,s2,-:,sL)  e  S,  set 

(2.1)  £(*)  =  {e€  E :  Sj  2  1  for  je  1(e)  and  Sj  —  0  for  je  L(e)\ 

so  that  E(s)  is  the  set  of  transitions  that  are  enabled  when  the  marking  of  the  spn  is  s. 
When  the  marking  of  the  spn  is  s  the  firing  of  an  enabled  transition  e€  E(s)  triggers  a 
marking  change  to  s'.  We  denote  by  p(s';s,e)  the  probability  that  the  new  marking  is  s' 
given  that  transition  e  fires  when  the  marking  is  s.  For  all  s  =  (jj,^.— .^). 
s'  »  (ijv^,-..,-*^)  €  S',  and  e€  E(s)  we  assume  that  p(s';s,e)  >  0  only  if 

(i)  Sj  -  1  *  sj  S,  sj  for  all  j€  1(e)  0  (D  -  J(e)), 

(ii)  Sj  —  1  £  s'j  £  Sj  +  1  for  all  je  1(e)  D  /(e), 

(iii)  Sj  s  s'j  £  sj  +  1  for  all  jeJ(e)  n  (D  -  7(e)),  and 

(iv)  s'j  m  Sj  for  all  je  (D  -  /(e)  -  7(e)). 

The  actual  enabled  transition  e  which  triggers  a  marking  change  when  the  marking  is  s 
depends  on  clocks  associated  with  the  enabled  transitions  and  the  speeds  at  which  these 


clocks  run.  Each  such  clock  records  the  remaining  time  until  the  transition  fires.  We 
denote  by  rsi  (  2  0)  the  deterministic  rate  at  which  the  clock  associated  with  transition  e( 

_ s  when  the  marking  is  s ;  for  each  s£  S,  rJ(  =  0  if  E(s).  We  assume  that  rsi  >  0  for 

some  et  e  E(s).  (Typically  in  applications,  all  speeds  rj(  are  equal  to  one.  There  are, 
however,  models  in  which  speeds  other  than  unity  as  well  as  state-dependent  speeds  are 
convenient.) 

For  s£  S  define  C(^)  to  be  the  set  of  possible  clock  readings  when  the  marking  is  s: 
C(j)  =  |(c1,...,Cjy):  Cj  2  0  and  c(  >  0  if  and  only  if  e,£  E(s ); 

(2.2)  c/"1  #  cjr-j1  for  »  j  with  CfCjrsirsj  >  o}. 

The  conditions  in  Equation  (2.2)  ensure  that  no  two  transitions  fire  simultaneously  as 
defined  below.  The  clock  with  reading  c,  is  said  to  be  active  when  the  marking  is  j  if 
transition  e(  is  enabled  (et  €  E(s)).  For  s€  S  and  ce  C(j),  let 

(2.3)  /*  *  /*(j,c)  =  min 

{i:et€E(j)} 

where  c/JJ1  is  taken  tobe  +*  when  rsi  =  0.  Also  set 

(2.4)  c]  *  c*(j,c)  «  c(.  -  t\s,c)rsi,  e,  €  E(s) 

and 

(2.5)  i*  -  i*(s,c)  —  i  such  that  £(i)  and  c*(i,c)  =  0. 

Beginning  with  marking  s  and  clock  vector  c,  t’(s,c)  is  the  time  to  the  next  transition 
firing  and  i*(r,c)  is  the  index  of  the  unique  firing  transition  e  =  e*(j,c)  =  cy 
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At  a  marking  change  from  s  to  s'  triggered  when  transition  e  fires,  new  clock  times 
are  generated  for  each  e  €N(s'\s,e  )  =  E(s)  —  ( E(s )  -  {**]).  The  distribution  function  of 
such  a  new  clock  time  is  denoted  by  F(-;r',e',j,e*)  and  we  assume  that  F(0;s\e  ,s,e*)  =  0. 
For  e  £  0(.s'-,s,e*)  =  E(s)  n  (£(*)  —  {?*}),  the  old  clock  reading  is  kept  after  e  fires.  For 
e'£(E(s)  -  {e*})  —  E(j'),  transition  e  (which  was  enabled  before  transition  e  fired)  is 
disabled. 

Next  consider  a  gssmc  z  0}  having  state  space 

2  =  U  (WxCW) 
ses 

and  representing  the  marking  ( Sn )  and  vector  ( Cn )  of  clock  readings  at  successive 
transition  firing  times.  (The  /th  coordinate  of  the  vector  Cn  is  denoted  by  Cnj.)  The 
transition  kernel  of  the  Markov  chain  {(Sn,Cn):/i  z  0}  is 

(2  6)  E((s,c)„4)  =  p(s\s,e  )  F(a|.;j',e|,j,e*)  1[0^  j(c*), 

e,eN(s)  e,EO(s ) 

where  N(s)  =  N(s;s,e*),  0(s )  =  0(j';r,e*),  and 

A  «  {/}  x  C(j'):  c]  z  a(.  for  e,.£  £(/)}. 

The  set  A  is  the  subset  of  Z  which  corresponds  to  the  spn  changing  marking  to  s'  with  the 
reading  cj  on  the  clock  associated  with  transition  e,£  E(s')  set  to  a  value  in  [0,a(].  (We 
suppose  that  the  clock  setting  distributions  are  such  that  F((s,c),S)  =  1.) 

Finally,  the  spn  is  a  piecewise  constant  continuous  time  process  constructed  from 
the  gssmc  {(S’n,Cn):n  z  0}  in  the  following  manner.  Set  =  0  and  denote  by  tn  the  time 
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of  the  nth  transition  firing,  n  *  0.  (We  assume  that 

/•{sup  =  +  »I(S,0,C0)}  -  1  a.s. 

n  i.  1 

for  all  initial  states  (S0,C0).)  Then  set 

(2.7)  X(t )  =  S' x(jy 
where 

(2.8)  #(/)  =  max  [n  i  0:  s  /}. 

The  process  {X(t):t  £  0}  defined  by  Equation  (2.7)  is  an  spn. 

Example  (2.9)  illustrates  the  use  of  the  spn  framework  for  formal  specification  of  a 
discrete-event  simulation.  For  ring  networks  with  #  ports,  reference  to  port  index  is 
to  be  interpreted  as  reference  to  index  j  -  1  (  mod  #)  +  1.  In  the  graphical 
representation  of  an  spn,  places  are  drawn  as  circles  and  transitions  as  bars.  Directed  arcs 
connect  transitions  to  output  places  and  input  places  to  transitions.  Tokens  are  drawn  as 
black  dots. 

(2.9)  EXAMPLE  (Token  ring  1).  Consider  a  unidirectional  ring  network  having  a  fixed 
number  of  ports,  labelled  1,2,...,#  in  the  direction  of  signal  propagation;  see  Figure  1.  At 
each  port  message  packets  arrive  according  to  a  random  process.  A  single  control  token 
(denoted  by  T  in  Figure  1)  circulates  around  the  ring  from  one  port  to  the  next.  The  time 
for  the  token  to  propagate  from  port  j  —  1  to  port  j  is  a  positive  constant,  R  j, 
j  *  1,2,...,#.  When  a  port  observes  the  token  and  there  is  a  packet  queued  for 
transmission  the  port  converts  the  token  to  a  connector  (C)  and  transmits  a  packet 
followed  by  the  token  pattern;  the  token  continues  to  propagate  if  there  is  no  packet 
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queued  for  transmission.  By  destroying  the  connector  prefix  the  port  removes  the 
transmitted  packet  when  it  returns  around  the  ring. 

Assume  that  the  time  for  port  j  to  transmit  a  packet  is  a  positive  random  variable, 
Lj,  with  finite  mean.  Also  assume  that  packets  arrive  at  individual  ports  randomly  and 
independently  of  each  other:  the  time  from  end  of  transmission  by  port  y  until  the  arrival 
of  the  next  packet  for  transmission  by  port  j  is  a  positive  random  variable,  Aj,  with  finite 
mean.  Note  that  there  is  at  most  one  packet  queued  for  transmission  at  any  time  at  any 
particular  port. 

Following  [4],  take  D  =  {1,2,. ..,4#}  and  E  =  [ex,e2,-^e3N\.  (See  Figure  2  for 
N  =  2.)  Set 

(2.10)  Z.(e3;_2)  =  UeZj_x)  =  L(e3j)  =  0, 

(2.11)  I(e3._2)  -  {4 j  -  2},  I{e3j_x)  =  {4 j  -  3,  4 j  -  1},  I(e3j)  =  {4(y  -  Dl, 
and 

(2.12)  /(e3y_2)  =  {4/  -  3},  J(e3j_x)  =  [Aj  -  2,  4 j],  J{e3j)  =  {Aj  -  1,  4yj, 

j  -  1,2 . N. 

The  transitions  have  the  following  interpretation:  e3;_2  =  "arrival  of  packet  for 
transmission  by  port  y,"  «3^_j  =  "end  of  transmission  by  port  y,"  and  e3  -  =  "observation 
of  token  by  port  y,"  j  =  1,2,. ...N.  The  interpretation  of  the  places  is  as  follows.  Let  12O 
and  s  *  (sx,s2,...,sAN)  €  S.  If  the  marking  of  the  spn  at  time  t  is  s  then  s4^_3  =  1  if  and 
only  if  at  time  t  port  y  is  transmitting  a  packet  or  there  is  a  packet  waiting  for 
transmission  by  port  y;  =  *  an<*  only  if  at  time  t  port  y  is  not  transmitting  a  packet 
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and  there  is  no  packet  waiting  for  transmission  by  port  j;  =  1  if  and  only  if  at  time  t 
port  j  is  transmitting  a  packet;  and  s 4 .  =  1  if  and  only  if  at  time  t  the  control  token  is 
propagating  to  port  j  +  1,  j  =  1,2,.. ..JV.  (Otherwise  Sj  =  0,  j  =  1,2,.. .,4 N.) 

Set 

S'  =  {(*lrr 2,...,j4JV):  sj  =  0  or  1  for  1  s  4W;  ^  +  s2  +  ...  +  s±N  =  N  +  1}. 

The  set,  S,  of  markings  is 

(2.13)  5  -  ((^^2,-,%)e^:  i4y_3  +  s4j—2  =  1  and  =  0  for  1  tf}. 

(It  follows  that  |  5 1  =  3N2n~*.  In  any  marking  there  are  exactly  #+l  tokens.  There  is 
at  most  one  token  in  each  place.  Each  of  the  disjoint  sets  of  places  {4/  —  3,4 j  —  2} 
contains  exactly  one  token  indicating  whether  or  not  port  j  has  a  packet  queued  for 
transmission.  The  set  of  places  D  -  {1,2,5,6,...,4jV  -  3,4N  -  2}  contains  exactly  one 
token  indicating  the  position  and  status  of  the  control  token.  There  can  never  be 
tokens  at  places  4 j  —  2  and  4y  —  1  simultaneously,  reflecting  the  fact  that  there  can  be 
no  arrival  of  a  packet  for  transmission  by  port  j  during  a  transmission  by  port  j.) 


The  new  marking  probabilities  are  as  follows.  If  e  =  e3^_2  =  "arrival  of  packet  for 
transmission  by  port  j"  then  /7($';j,e)  =  1  when 

s  =■  (Ji,...^4^_jj,0,l,0,J4y,i4y+1,...,r4^)  £  51  and  s  —  (rj,...,r4^_j^,l ,0,0,r4y,r4y+1,...,r4^). 

If  e  —  e3 j_l  =  "end  of  transmission  by  port  j,"  then  p(s'\s,e )  =  1  when 

s  **  (^••••44y.jjil>0,l,0^4y+i,.-  rf4^)£  ^  and  s  —  (jj,...,j^^_j^,0,l,0,l,J4y+j,...,i4^y). 

If  e  -  e3y.  *  "observation  of  token  by  port  j”  then  p(s'\s,e)  =  1  when 
J*  (rji...rs4^j)_jil,li0,0,0,i4y+j,...,...^4^)£  S  and  s  =  (•rj,-..,J^^_j^_j,0,l,0,l,0,j^+ 


k  ■  J  ■  V 
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and 

T/0  =  inf  {{•„  >  S\_,(/):  *(£„)  e  ^(^_i)  e  *iL  J  *  *• 

Note  that  since  the  spn  is  symmetric  the  sets  A^,  A l2,  B^t  and  satisfy  the  conditions 
which  ensure  that  the  start  and  termination  times  for  the  passage  times 
Pj  =  Tj(f)  —  S\_j  (/)  strictly  alternate.  Denote  the  successive  passage  times 
p\,p\,...,P^,p\,p\t...  enumerated  in  termination  order  by  {P'jij  i  1}.  Set  7^  =  0  and  let  7 ^ 
be  the  termination  time  for  Pj,  j  k  1. 

Proposition  (4.1)  gives  conditions  which  ensure  that  {(Z( Jj|),/’)|+1):n  i  0{  is  a 
regenerative  process  in  discrete  time  and  that  the  expected  time  between  regeneration 
points  is  finite.  The  regenerative  structure  guarantees  (Miller  [13])  that  P'n^P  as  n-*». 
The  goal  of  the  simulation  is  the  estimation  of  r(f)  *  £[/(P)},  where  /  is  a  real-valued 
(measurable)  function  and  P  is  the  limiting  passage  time. 

We  postulate  the  existence  of  a  transition  e*  and  markings  *0e^i  an<*  Joe^2  suc^ 
that  a  passage  time  terminates  when  transition  e  fires  and  the  marking  changes  from 
to  jq-  Io  addition,  we  assume  that  no  passage  times  are  underway  when  the  marking  of 
the  spn  is  Jq.  Formally,  denote  by  L(t)  the  last  marking  of  the  spn  before  jumping  to  X(t) 
and  set 


V{t)  -  (L(r)  ,*(/))• 

Denote  by  G  the  state  space  of  {V(t):t  z  0}.  Set  A1  =  A^xA^  and  B1  =  B^xB^, 

//j  =  {j  €  5:  (s,s  )  €  Bl  -  A1  for  some  jes} 


/  =  1,2,...,//.  Now  set 
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4.  PASSAGE  TIMES  IN  STOCHASTIC  PETRI  NETS 

Formal  specification  of  passage  times  in  a  symmetric  spn  is  by  means  of  four  subsets 
(Aj,  A2,  B},  and  B2 )  of  the  marking  set,  S.  The  sets  Aj,  A2,  B j,  and  B2  in  effect 
determine  when  to  start  and  stop  the  clock  measuring  i  particular  passage  time;  cf. 
Iglehart  and  Shedier  [7]. 

Denoting  the  jump  times  of  the  process  { X(t):t  £  0}  by  ^  0},  for  k,n  i  1  we 

require  that  the  sets  Alt  A2,  B  j,  and  B2  satisfy: 

if  X{iH_l)€Av  *tt„)£A2,  X($n_l+k)£Av  and  *(fn+Jfc)eA2 
then  X($n_l+m)e  and  X(Sn+m)E  B2  for  some  0  <  m  s  k\ 

and 

if  A(fn_,)  e  B  j ,  *(£„)  e  B2,  +k)eBv  and  *(fn+*)  e  B2 

then  X(Sn_l+m)eAx  and  X(£n+m)EA2  for  some  0  £  m  <  k. 

These  conditions  ensure  that  the  start  and  termination  times  for  the  specified  passage 
time  strictly  alternate. 

In  terms  of  the  sets  Av  A2,  Bv  and  B2 ,  define  Alx  =  {<#>5(5):  sEA^}, 
A2  *  s£A2],  sEB j},  and  B2  =  {^(s):  sEB2],  l  =  1,2,. ..,N.  (Recall 

that  for  sES,  ^(j)  =  <f>s(s )  and  <t>ls(s)  =  (*)) .  /  =  1,2 Also  recall  that  for 

e£  E,  $£(«)  -  4>£(e)  and  <t>E(e)  =  <f>E(.<f>lEX(e)).)  Then  define  two  sequences  of  random 
times  {Sj(l)'J  *  0}  and  {Tj([):j  2  1}:  S\_ t ( /)  is  the  start  time  for  the  /th  passage  time 
(corresponding  to  the  sets  A Al2,  flj,  and  Bl2 )  and  Tj(l)  is  the  termination  time  of  this 
/th  passage  time.  Set 

S0(l )  =  0, 

Sjd)  =  inf  {{„  2  TjV):  X^EA1^  *(fn_j )EA\},j  *  1 
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Standard  arguments  establish  a  ratio  formula  for  r(J). 


(3.15)  PROPOSITION.  Provided  that  E[rx}  <  «  and  E{  |/U)  | }  <  *, 


rif)  =  E[f{X)\ 


With  these  results  Equations  (2.22)  and  (2.23)  provide  point  estimates  and  confidence 
intervals  for  r(/). 


(3.16)  EXAMPLE.  In  the  token  ring  model  of  Example  (2.18),  take  <f>(J )  =  j  +  1, 
j  =  1,2,...,N.  Set 

S'  —  { seS :  s =  1  for  some  j,  j  =  1,2,...,^} 

and  consider  the  function  /  defined  by 

/W  -  lj5'}(*) 

for  s€  S.  According  to  this  definition,  r{f)  is  the  steady  state  throughput  of  the  token 
ring.  Note  that  the  function  /  satisfies  Equation  (3.1)  since  (for  each  0  _ i  =  1  if 

and  only  if  =  1,7=  1,2,..., N.  Arguments  given  in  [4]  show  that  the  conditions 

of  Proposition  (2.19)  hold  provided  that  the  packet  interarrival  time  random  variables, 
Aj,  have  new  better  than  used  distributions  and  satisfy  a  positivity  condition: 
P[Aj  s  Rn\  >  0,  y  =  1,2 
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-  P,^  =  n  +  1,  /(*(£„))  =  /(*„),  SH  i  zn, 

/w»rl))  =  s  =  ♦*<*>} 

for  all  i  0,  S,  and  nil.  Applying  Equation  (3.14)  and  then  Equation 

(3.13),  it  follows  that  for  all  Zj,...^  i  0  and  jclt...,xne  5: 

M«*+i  -  »  +  WTlJ)  =  /(V.  fy  +„  *  v 

•^■*^yk+n-l^  “  Xn -P*  ^yk+n-l  S  Z/i-r"'’^y*+l  5  zl^ 

N 

-  2  Wi  -  »  +  1.  W)  -  A*.).  fB  *  v 

/-i 

~  A*,_ j).  f„_j  s  *,_!»— t£i  *  zx}  /*i{*(fyfc)  -  <#>$(*())} 

N 

-2  ^it*i  -  »  +  1.  AW*))  =  A*„).  f„  *  zn, 

/-i 

~  A-*n_j).  *  . . . ?1  s  *l)  ^i{-^(fyk)  =  4>5(jo^} 

-*,{«!-  n  +  l.  AX(S„))  =  /(*„),  rn  *  Z„,  A^f^j))  -  A*^).  f„_!  £  Vi . h  *  zi* 

so  that  (using  Equation  (3.10))  the  pairs  of  random  variables  {(Yk(f),rk):k  i  1}  are 
identically  distributed.  □ 
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'*♦1  '  r»,.,  -  Th  and  the  finite  dimensional  distributions  of  {/(AT(/)):/  2  V-  Observe 
that  the  joint  distribution  of  *(7^)  =  $$(sq)  and  the  clocks  set  or  reset  at  time 
depends  on  the  past  history  of  the  spn  only  through  the  new  marking  the  previous 

marking  $$( sQ ),  and  the  trigger  transition  ^(e*);  this  implies  that  the  cycle  length  r;fc+1 
and  {/(AX/)):/  2  are  independent  of  l(Yj(f),Tj):j  £  A}.  It  follows  that  the  pairs  of 
random  variables  {(Yk(f),rk):k  2  1}  are  mutually  independent. 


Next  observe  that 

r*-1 

(3.10)  W-  £  /WW)  trm+1  -  U. 

where  {  *  7^,  k  2  0.  Set  Sk  =  yk  -  yk_v  k  i  1.  It  is  sufficient  to  show  that  for  all 

2  0,  xv...jcne  S,  and  nil: 

«  +  l,  m {„))  -  /<*„).  ?„  s  v -  A*„-1).ril_1  s  z„_, . t,  ^  *,) 


(3.11)  -  -  »  +  1.  A*fr,.„)>  -  A*„).  fTt„  *  I„  AWr, *  A*..,). 


^y4+fl-i  s  Irk+rt-r  '  ^7fc+i  s  zi 

Ail.  Here  denotes  the  conditional  probability  associated  with  starting  the  spn 

with  marking  ^(jq)  and  all  active  clocks  reset  at  time  /  =  0  according  to  the 
distributions 

(3.12)  P{CQ' «•  5  c)  »  F(c;$^(jQ),el-,^y(jg),££(e  )) 

for  c  i  0,  eie  £($^(jq)).  Recall  that  CQ.  is  the  clock  reading  associated  with  transition  e, 
at  time  0.  Denoting  by  Pt[-}  the  corresponding  conditional  probability  when  the  initial 


P{X{ 7i)  -  uroi.j)  -  =  ♦iUj)}, 


/  —  1,2,.. .,N.  The  result  follows  from  Equation  (3.3).  □ 

Carry  out  the  simulation  of  {A'(r) :/  z  0}  in  random  blocks  defined  by  the  successive 
random  times  {7^:*  z  0},  defined  by 

(3.8)  7„  -  inf {7^  >  Tb  :  *(7^,)  »  <t>‘s(s'0)  for  some  /,  /  =  1,2,...,*}, 

ill;  P0  *  0  and  7^  =  0.  (Note  that  the  random  times  {7^:/c  *  0}  do  no/  form  a 
sequence  of  regeneration  points  for  the  process  (Jf(/):/  z  0}.) 


ill. 


(3.9)  PROPOSITION.  The  sequence  of  pairs  of  random  variables  {(.Yk(f),rk):k  2:  1}  are 
independent  and  identically  distributed. 


Proof:  The  sequence  {Pk:k  i  lj  are  indices  of  the  successive  stopping  times  { Tn:n  *  1}  at 
which  transition  fires  and  the  marking  changes  from  $$(/)  to  <^( sQ )  for  some 

s  eS  and  some  /,  /  =  1,2,...,*.  Thus,  by  the  definition  of  a  symmetric  spn,  each  of  the 
clocks  running  at  time  Tp +  was  set  or  can  be  viewed  as  having  been  probabilistically 
reset  at  time  TK  Therefore  [X(t)u  z  V  determines  the  distribution  of 


(3.5)  P{X(SH)  =  4(^),  en  *  =  4^),  e\  -  4^).  *(f0)  =  4<;c0)} 

=  P^X{^  —  $s(xn),  en  —  ^£(<I|()*-I^(fi)  =  4^*l)*  el  ~  4(e»',)»  X({q)  —  4^x0^}’ 

I  =  l,2,...,Af.  Denote  by  { yln:n  2  0}  the  indices  of  the  successive  times  tf„:n  z  0}  at  which 
transition  4>E(e)  fires  when  the  marking  is  4>ls(s*)  for  some  s*eS*.  Equation  (3.5) 


•  •  _• 


implies  that  for  all  xx,...jcn£  S  and  sx,...tsne  S  : 


(3.6)  W*(ryo  -  4<*„>’  e\  -  wr<_!)  =  4  ■<£.•••.*»/)  =  4<*i>» 


e\  =  <f>lE(e),  =  4(4)} 

=  P{XU  ,)  =  4(jr„),  e\  =  4(e*),  Wj-J  =  4 =  4<*1>* 

*1  “  4^  )»  *  4^l)}. 

/  =  1,2,...,M  Using  the  definition  of  {7^:n  Z  0},  Equation  (3.6)  implies  that  for  all 
xlf...,xn€  S: 

(3.7)  Ff^rJ)  =  4u„),  ...,*( r})  -  4(Xi)}  *  F{^(4)  =  4(x„),  ...^(4)  =  4(xj)} 
/  =  1,2,.. .,Af.  Applying  Equation  (3.7)  it  follows  that  for  all  x,  ,  €  S: 

*  /l“  1 

F{^(ri)  «  40q)  MKr^)  =  ^(v!).-^^  =  4<*i>} 

p{x{t\)  «  4(^q),  x{t\_x)  -  4(xn_j),...^(rJ)  =  4ut)} 


'{*<4-1)  *  4<*„-i>"  ^ri>  -  4<*i)} 

p{x{  4) = 4(5q),  = 4^n-i)-^4)  -  4<*i>} 


F{^(4_i)  =  ♦i(V_l)..  .^f(4)  -  4(^)} 
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3.  STEADY  STATE  ESTIMATION  FOR  SYMMETRIC  STOCHASTIC  PETRI  NETS 

In  this  section  we  consider  estimation  of  rif)  under  the  assumption  that  the 
function  /  is  symmetric  in  the  sense  that 

(3.1)  /W  -  /t4W) 

for  all  s£S  and  all  l  —  1,2,...,N.  Symmetry  of  the  underlying  spn  implies  that 
regenerative  cycles  defined  by  the  times  at  which  the  transition  ^(e*)  fires  and  the 
marking  changes  to  <^(5q)  can  be  decomposed  into  independent,  nonidentically 
distributed  blocks.  These  blocks  are  defined  by  the  successive  times  Tn  at  which 
transition  <t>lE(e*)  fires  and  the  marking  changes  from  <#>^(/)  to  $ls(s0)  for  some  /  e  S *  and 
some  /,  /  —  1,2 Estimates  for  K/)  can  be  based  on  observation  of  these  blocks. 
Proposition  (3.2)  provides  conditions  which  ensure  that  (for  each  0  transition  4>lE(e) 
fires  and  the  marking  changes  to  <t>ls(s'Q)  infinitely  often  with  probability  one. 

Denote  by  { Tn:n  2  1}  the  times  ,T\,—  in  increasing  order. 

(3.2)  PROPOSITION.  Suppose  there  exists  fi  >  0  such  that 

(3.3)  />{Jf(rJ)  -  4(^)|jf(r^1),...^(rJ)}  *  s  a.s.. 

>  Then  P[X{Tn)  -  «/5( sQ )  i.o.}  =  1  for  all  /  =  1,2,...,N. 

Proof:  By  Lemma  4  of  [6]  it  suffices  to  show  that 

>  (3.4)  /»{A-(Ti)  =  4(io)U(7Vn_1),...^(r,1)}  2  S  a.s. 


Let  e*  denote  the  transition  that  fires  at  time  £n,  n  2  0.  The  definition  of  a  symmetric 
spn  implies  that  for  all  S'  and  e,  e  E: 


15 


consistent  point  estimate 

(2.22) 


Arin) 


Y{n) 

T  («) 


and  asymptotic  100(1  —  2y)%  confidence  interval 


(2.23) 


*l_y  *(") 

r(n)  nV1 ' 


r(n)  + 


rl_7  *(«)  " 

?(«)  nV1  _ 


for  r(/).  In  Equation  (2.22) 

n«)  =  «_1  £ 

m  — 1 


and 

=  n_1  £  Tm> 

m-1 


where  Tm  is  the  length  of  the  /nth  cycle  and  Ym(f)  is  the  integral  of  over  the  mth 

cycle.  The  quantity  ,s2(n)  is  a  strongly  consistent  point  estimate  for 

a2(J)  =  var  (Tj(/)  -  K/)tx) 

and  =  «h-1(l  —  y),  where  is  the  distribution  function  of  a  standardized  normal 
random  variable,  N(0,1).  Confidence  intervals  are  based  on  the  central  limit  theorem 
(c.l.t.) 


n1/2  }r(n)  -  r(J)} 

Oif)/ E[  Tj} 


mi) 


as  n-* oo. 


Also  suppose  that  there  exists  seS  such  that  for  all  s  £  S  , 

(i)  the  set  O(sqSS)  -  n  (£(/)  -  {/})  =  0, 

(ii)  the  set  N(sQ\s*,e*)  =  E(s^)  -  ( E(s *)  -  {**})  -  N^;*,/),  and 

(iii)  the  clock  setting  distribution  F(;s0,e,s*,e*)-F(;s0,e\s,e*)  for  all  e  e  N(s0;s,e*). 


Then  {*(/):/  a  0}  is  a  regenerative  process  in  continuous  time.  Moreover,  if 

(2.21)  e{t\+1  -  rj}  s  c  <  - 

for  all  n  i  0  then  the  expected  time  between  regeneration  points  is  finite. 

Equation  (2.20)  implies  that  transition  ^(e*)  triggers  a  marking  change  to  ^(jq) 
infinitely  often  with  probability  one.  Furthermore,  at  such  a  time  T\,  the  only  clocks 
that  are  active  have  just  been  set  since  0($£(5q);<^(/), $£•(/))  —  0  for  all  /  £  S* .  The 
joint  distribution  oi  X(t\)  and  the  clocks  set  at  time  7*  depends  on  the  past  history  of 
{X(/):t  2  0}  only  through  <^(jq),  the  previous  marking  £$(/),  and  the  trigger  transition 
$£(<?*).  Since  the  new  transitions  and  clock  setting  distributions  are  the  same  for  all  /, 
the  process  i  0}  probabilistically  restarts  whenever  {X(T*):n  2  1}  hits  ^(jq). 

Note  that  the  result  of  Proposition  (2.19)  also  holds  if  condition  (i)  is  replaced  by: 

(i')  O(jo;jo,e*)f»0  and  for  any  e£  0(jQ;jQ,e*)  the  clock  setting  distribution  F(  ;s,e\s,e)  is 
exponential  with  mean  which  is  independent  of  s,  s',  and  e.  (Assumption  (i1)  ensures  that 
no  matter  when  the  clock  for  transition  e'e  0(jq;jo,/)  was  set,  the  remaining  time  until 
transition  e  triggers  a  marking  change  is  exponentially  distributed  with  the  same  mean.) 

Under  the  conditions  of  Proposition  (2.19),  the  basic  limit  theorem  for  regenerative 
processes  asserts  that  X(t)+X  as  /-»*.  The  goal  of  the  simulation  is  the  estimation  of 
r{f)  —  E\J\X)},  where  /  is  a  real-valued  (measurable)  function  having  domain  S.  From  n 
cycles  the  standard  regenerative  method  (Crane  and  Iglehart  [2])  provides  the  strongly 


LvL2t—,Lff  are  identically  distributed.  Under  these  assumptions  the  spn  is  symmetric. 
For  example,  take  $(/)  =  j  +  1,  j  =  1,2,...,#.  Then 

*0(4 O'  -  1)  +  *)  -  4 j  +  *, 


*£**3(/-l)+P  “  e3j+k* 


and 


Proposition  (2.19)  (cf.  Proposition  (4.7)  of  [4])  gives  a  set  of  conditions  on  the 
building  blocks  of  an  spn  which  ensure  that  the  process  is  regenerative  and  that  the 
expected  time  between  regeneration  points  is  finite.  Set  =  <t>s(s)  and 

$$(*)  =  ^5(^5  *($)) 

for  s€  S  and  /  —  2 ,...,#.  Similarly,  set  $^(e)  =  $£(e)  and 

♦ito  *  ♦*(*Jr1(*)) 

for  «e£.  Recall  that  is  the  nth  transition  firing  time,  n  i  0.  Let  {7j,:n  i  0}  be  an 
increasing  sequence  of  stopping  times  that  are  finite  (7^  <  *>  a.s.)  transition  firing  times 
such  that  for  some  e  e  E  and  S*sS:  Tq  =  0  and 

ij,  =  inf{/  >  7^_j:  at  time  r  transition  0^(e")  fires  and  the  marking  is  <^(/)  for  some  j*e  S'* 
nil  and  /  =»  1,2,...,#. 

(2.19)  PROPOSITION.  Suppose  that  there  exists  Jq£  S  and  S  >  0  such  that 


(2.20) 


P{X(T\)  -  4^0)  U(rJ_1),...^f(r{)}  *  5  a.s. . 


12 


i 

i 

i. 

s 


Let  $  be  a  cyclic  permutation  of  the  set  {1,2,...,N}.  In  terms  of  this  permutation 
define  a  mapping,  <t>D,  of  D  onto  D : 

(2.14)  ♦*((/  -  DLj  +  *)  =  (*0)  -  D^i  +  *, 

y  *•  1,2,.. .Jf  and  k  -  l,2,...,Lj.  Similarly,  define  a  mapping,  $£,  of  E  onto  E : 

♦^(eo— l W,  +*)  “  %U)-lW,+*‘ 

y  =  1,2,.. .,N  and  k  =  1,2,..., A/,.  Also  define  a  mapping,  of  S'  onto  5: 

(2.16)  ^5(^lrf2,...^L)  *  (^(1)^X>(2).-^Z>(L))- 

For  O'  c  O,  we  write  <t>D(D)  =  {^(i):  ieO'}  and  for  E'  c  £  we  write 
$e(£  )  m  {$£(«)”>  *€  E  }. 

(2.17)  DEFINITION.  An  spn  {X(/):/  *  0)  is  said  to  be  symmetric  if  there  exists  a  cyclic 
permutation,  of  the  set  (1,2,.. „N]  such  that: 

(i)  <t>D{Ue))  =  U4>Eie)),  <f>DU(e))  =  /(*£(e)),  and  ^(/(e))  =  J(<f>E(e ))  for  all  e€  E, 

(ii)  p(s;s,e)  *  p(<frs(s);4>s(s),<t>E(e))  for  all  e€  E  and  s,s  e  S,  and 

(iii)  F(;s,e,s,e)  =  F(;$s(i'),<fr£(e'),<f»5(*),4>£(e))  for  all  e  e  N(j';s,e),  e€  £,  and  s,s  €  S. 

Condition  (i)  ensures  that  the  induced  mappings  <l>D  and  <pE  preserve  the  sets  of  normal 
input  places,  inhibitor  input  places,  and  output  places  and  thus  (using  Equation  (2.1)) 
$e(E(s))  —  E(<t>s(s))  for  all  se  S.  Conditions  (ii)  and  (iii)  ensure  that  the  mappings  <ps 
and  $£  preserve  the  new  marking  probabilities  and  the  clock  setting  distributions. 

(2.18)  EXAMPLE.  In  the  token  ring  model  of  Example  (2.9),  let  N  be  the  number  of 

ports  so  that  Lx  -  4  and  Mx  -  3.  Suppose  that  (i)  *  R2  »  ...  -  RN\  (ii)  the  random 

variables  A^2,...,AN  are  identically  distributed;  and  (iii)  the  random  variables 


r 


k 


t 
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and  when 

s  -  (s1.-.J4(/_1)_1,l,0,l,0,0,i4/+1 . sAN)£S  and  s  =  (s1,...^4^_1)_1,0, 0, 1,0,1, j4y+1,...,j4Ar). 

All  other  new  marking  probabilities  p(s\s,e)  are  equal  to  zero. 

Note  that  when  transition  e ^  fires,  a  token  is  removed  from  place  4 (J  —  1)  and  a 
token  is  deposited  either  in  place  4/  -  1  or  in  place  4 j,  depending  upon  whether 
(i4^_3rf4^_2)  equals  (1,0)  or  (0,1).  All  other  transitions  are  input-deterministic  (in  that 
exactly  one  token  is  removed  from  each  input  place  when  the  transition  fires)  and 
output-deterministic  (exactly  one  token  is  deposited  in  each  output  place  when  the 
transition  fires). 

The  distribution  functions  of  new  clock  times  for  transitions  e  e  N(s;s,e')  are  as 
follows.  If  e  =  e3y_2  *  "arrival  of  packet  for  transmission  by  port  j,"  then  the 
distribution  function  F(x;s  ,e  ,s,e")  =  P{Aj  £  x}.  If  e  =  e3^_1  =  "end  of  transmission  by 
port  j,"  then  the  distribution  function  /‘(x;j',e',j,e*)  =  P{Lj  £  x}.  If  e  =  e^  -  — 
"observation  of  token  by  port  j,"  then  the  distribution  function 
£(x;j' =  lf;?  ,„)(*)• 

We  now  define  a  symmetric  spn.  Informally,  an  spn  is  symmetric  if  there  are 
mappings  of  places  onto  places,  markings  onto  markings,  and  transitions  onto  transitions 
which  preserve  the  sets  E(s)  of  enabled  transitions,  the  new  marking  probabilities 
p(s';s,e),  the  sets  N(s;s,em)  of  new  transitions,  and  the  clock  setting  distributions 
£(•;*',*',*,**).  Let  {X(t):t  a  0}  be  an  spn  with  finite  marking  set,  S,  and  transition  set,  E. 
Throughout  this  section  we  assume  that  D  =  {1,2,...,L}  is  the  index  set  of  places  and 
£  -  is  the  transition  set,  where  L  =»  and  M  =  M^N  for  some  N  £  2. 

(We  assume  that  all  clock  setting  distributions  have  finite  mean.) 


so  that  h\  is  the  set  of  all  possible  markings  when  a  passage  time  Pln  terminates.  Also  set 


=  |j'  £  S:  for  all  se  S,(s,s)  £  G  -  (B1  U  A1  )  and  (*,*'), 

A1  i 

2)  for  some  (s^s^eB1,  (s^s^eA1] 

so  that  H2  is  the  set  of  all  possible  markings  when  a  passage  time  7^  is  not  underway. 

Ai 

(For  (j,5'),(5,j) 6  G  we  write  (s,s')~(j,j)  if  there  exists  a  finite  sequence  e0,s\,e\,s2,...jn,en 
of  transitions  and  markings  such  that 

and  (s'ts'1),(s^t,s),(sj,sj+1)  ^Al,  j  =  1,2,.. —  1.)  We  assume  that 

b2  =  b2  n  (#}  u  h\)  n  ...  n  (7/f  u  h%)  *  0 

and  that  sQe  B2. 

As  in  Section  2,  let  {7 Jn:n  i  0}  be  an  increasing  sequence  of  stopping  times  that  are 
finite  (7^  <  *>  a.s.)  transition  firing  times  such  that  for  some  e  €  E  and  S*zS:  7^  *  0  and 

7^  =  inf  {/  >  7^_j:  at  time  t  transition  4>lE(e*)  fires  and  the  marking  is  for  some  s'  e 

nil  and  /  =*  1,2,...,W. 

(4.1)  PROPOSITION.  Suppose  that  there  exists  e  e  E,  sQeBv  and  jgeflj  SUC^  t^lat 
p(sQ-fs0,e')  >  0  and  either  (i)  0(jq;jo,?*)  =  £(sq)  n  (£(*0)  -  {e*})  =  0  or  (ii) 

0(sQ;sQ,e)  0  and  for  any  e  £  £?(jQ;j0,e’)  the  clock  setting  distribution  F(-;s\e\s,e )  is 
exponential,  independent  of  s,  s,  and  e.  Set  Vq  =  (<#eI (^0) ,<#»^(^q))  and  suppose  there  exists 
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8  >  0  such  that 

(4.2)  P{V(Tln)  -  vj|  nr^) . nrj)}  a  8  a.s.. 

Then  {(XCfJtP'^-^'.n  2  0}  is  a  regenerative  process  in  discrete  time.  Moreover,  if 

-  Tl]  *  c  <  - 

for  all  n  2  1  then  the  expected  time  between  regeneration  points  is  finite. 

Proof:  Since  Tln  <  «  a.s.  and  P{V(Tln)  =  |  FfTj^) . V(T\)}  2  8  >  0,  Lemma  4  of  [6] 

ensures  that  transition  <pE(e  )  fires  and  the  marking  of  the  spn  changes  from  to 

infinitely  often  with  probability  one:  P{V(TX)  =  Vq  i.o.}  =  1.  Denote  by 
{/?*:*  2  1}  the  indices  of  the  successive  passage  times  { P'n:n  2  1}  which  terminate  when 
transition  )  fires  and  the  marking  changes  from  <*>^(s0)  to  Let  7q  *  0q  -  0. 

We  must  show  that 

(i)  iP\'-k  2  0}  is  a  renewal  process  in  discrete  time 


and  that  for  any  ix  <  t2  <  •••  <  im  (m  2  1)  and  *20 


(ii)  and  ha.e 


the  same  distribution,  and  Wf.t  .  ,f  is 

Pk+I  J  Pk+1 1+1  Pk+,m  Pk+,m  +  l 

independent  of  {(A(7^),/’J|+1):0  s  n  <  fi\}. 


At  time  7^i,  a  passage  time  has  just  terminated  with  no  other  passage  times  underway. 

Now  observe  that  each  of  the  clocks  running  at  time  Tf„ i  +  was  set  or  can  be  viewed  as 

Pk 

having  been  probabilistically  reset  at  time  7^i.  (Assumption  (ii)  ensures  that  no  matter 
when  the  clock  for  transition  e  6  0($j(jQ);$].(,so),^(e  ))  was  set,  the  remaining  time 
until  transition  e  fires  is  exponentially  distributed  with  the  same  parameter.)  Therefore 


{*(/):/  2  7^}  determines  the  finite  dimensional  distributions  of  ^(7^i+/),  ^i+f+1  for 
i  2  0  and  the  distribution  of  /}£+1  -  The  joint  distribution  of  *(7^1 )  and  the  clocks 
set  or  reset  at  time  7^i  depends  on  the  past  history  of  the  spn  only  through  the 

previous  marking  ^(r0),  and  the  trigger  transition  <t>E(e  ).  This  distribution  is  the  same 
for  all  and  therefore  (i)  and  (ii)  hold. 

Proposition  (4.3)  of  [4]  implies  that  2  0}  is  a  regenerative  process  in 

continuous  time  and  , ,  —  7^_i }  <  «.  It  follows,  since  the  state  space  of  the  spn  is 

Pk+1  Pk 

finite  and  the  clock  setting  distributions  have  finite  mean,  that  Eip\+\  —  <  ».  □ 

Proposition  (4.3)  provides  sufficient  conditions  which  ensure  that  Equation  (4.2) 
holds.  We  postulate  the  existence  of  a  distinguished  random  time  Tf  in  the  interval 
[7,*_1.7’*)  and  a  set  [e^:k€  Jf(v+)}  of  distinguished  transitions  determined  by  the 
marking,  v+,  at  time  T%.  We  make  the  following  sample  path  assumption :  V(TXn)  =  vj  when 
each  of  the  distinguished  transitions  occurs  prior  to  some  time  +  ^^(v*). 
Proposition  (4.3)  asserts  that  the  geometric  trials  recurrence  criterion  (Equation  (4.6))  is 
satisfied  if  the  clock  setting  distributions  associated  with  the  distinguished  transitions  are 
"new  better  than  used"  (NBU)  and’satisfy  a  "positivity"  condition  (condition  (iii))  which 
guarantees  the  existence  of  5  >  0  as  in  Equation  (4.6).  (A  positive  random  variable  A  is 
NBU  if 

P{A  >  x  +  y  \A  >  y)  z  P{A  >  x\ 

for  all  xy  2  0.  Note  that  every  increasing  failure  rate  (IFR)  distribution  is  NBU.  Also, 
if  A  and  B  are  independent  random  variables  with  NBU  distributions,  then  the 
distributions  of  A  +  B,  min  ( A.B ),  and  max  (A,B)  are  NBU.) 


Recall  that  G  is  the  state  space  of  the  process  {V(t):t  a  0}.  Let  {T%:n  t  0}  be  a 
sequence  of  transition  firing  times  and  denote  the  state  space  of  {V(T^):n  z.  0}  by  G+. 
Set  —  {(S/tC/):0  s  /  <  M^)}.  where  N(-)  is  given  by  Equation  (2.8).  Let 

e(1),e(2) . «(m)£  E  and  for  v+  -  (/+^e+)£  (7+,  set  £(v+)  -  £(x+ )  and 

K(v+)  -  {*:  e(W££(v+)}. 

When  y(T^)  —  v+,  for  k £  /f(v+)  we  denote  by  SnJc(v+)  the  latest  time  less  than  or  equal 
to  at  which  the  clock  associated  with  transition  e ^  was  set  and  by  A^^v*)  the 
setting  on  the  clock  at  time 

(4.3)  PROPOSITION.  Let  £  E  and  let  {TJ^rniO}  be  a  sequence  of 

transition  firing  times.  For  v+  £  G+,  let  {jRn^(v+):fc  £  /f(v+)},  be  identically  distributed 
collections  of  random  variables,  independent  of  [AnJ((V{T^)):ke  K(V(T^))}  and  ^jf(r^). 
Assume  that: 

(i)  a.s.  and  for  G  and  v+£  (7+, 

(4.4)  p{v<t\)  -  vj,  no  -  v+.  nr'_,)  =  Vl . nr!)  =  v0} 

2  '■{^(v+)  +  A'j'iv+l  s  7 *  +  RnJt(v+).  ((«.*); 

no  -  -+.  nr'.,)  -  Vl . nr])  -  v0}; 

(ii)  for  ail  e(/c)  the  clock  setting  distribution  F(;j',e(*\j,e)  =  F(  ;e(A:))  and  is  NBU;  and 


(iii)  there  exists  8  >  0  such  that  for  v+  £  G+ 


(4.5)  fi(v+)  -  H4k(v+)  S  RnJl(v+),  k  £  AXv+)  J  i  8, 

where  the  random  variable  A^v+)  has  distribution  F(  ;e^)  and  L4{v+):/£  JC(  v+)}  are 
mutually  independent  and  independent  of  ,(v+):/£  AT(v+)}. 

"v 
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Then 

(4.6)  P{V{Txn)  =  vj  |  nrJ_1),...,K(r;)}  *  8  a.s 

so  that  P[V(T\)  =  vj  i.o.}  =  1. 

Proposition  (4.3)  follows  directly  from  Proposition  (2.16)  of  [3]  since  the  process 
2  OJ  is  a  generalized  semi-Markov  process  with  state  space,  G,  and  event  set,  E. 


(4.8)  EXAMPLE.  In  the  token  ring  model  of  Example  (2.17),  take  $(y)  =  j  +  1, 

t 

j  —  1,2,..., AT.  Set  s0  =  (1,0, 0,0,... ,1,0, 0,0, 1,0, 0,1)  and  s'0  =  (1,0,1, 0,1, 0,0,0,...,1, 0,0,0). 
Take  e  =  e3  and 


5  —  {(ij^i— ^jy)  £  5:  ^  =  1 1 


(where  S  is  given  in  Equation  (2.13))  so  that  t\  is  the  nth  time  at  which  port  2  observes 
the  token,  nil.  (Note  that  =  $j(jq)  if  there  is  a  packet  queued  for  transmission 

at  each  of  the  other  ports  and  port  2  starts  transmission  of  a  packet  at  time  T} J.  The  spn 
{X(i):t  2  0}  changes  marking  to  ^(jg)  when  transition  )  fires  and  the  current 

marking  is  <^(*0).)  Observe  that  T*  <  «  a.s.  since 


E{Txn  -  rJ_j)  E\LX }  <  oo 


for  all  n  1.  Take  =  e3y_2  (m  »  N).  Let  be  the  first  time  after  at  which 
the  control  token  leaves  port  1  (transition  )  becomes  enabled).  Take  ^„^(*'+)  =  P\ 
for  all  v+e(7+.  Since  the  spn  has  marking  <p\(s'q)  at  time  t\  if  each  transition  e3^_2 
enabled  at  time  fires  before  <J>£(e3)  fires,  condition  (i)  of  Proposition  (4.3)  is  satisfied. 
Assume  that  for  j  =  1,2 (i)  the  distribution  of  Aj  is  NBU  and  (ii) 


* 


Sj  —  s  /?  j}  >  0  so  that 

N 

«(v+)  -  n  *j  *  n  sj  - 6  >  °* 

j€K(*+)  ;"1 

Then  i*W7j)  -  vj  i.o.}  -  1. 

The  definition  of  a  symmetric  spn  implies  that,  for  the  process 
{(A'(7^)^P)|+1):n  2  0},  regenerative  cycles  defined  by  the  times  at  which  the  transition 
$£(e  )  fires  and  the  marking  changes  from  <p^(sQ)  to  <f>^(s0)  can  be  decomposed  into 
independent,  nonidentically  distributed  blocks.  These  blocks  are  defined  by  the 
successive  times  Tn  at  which  transition  ^(e*)  fires  and  the  marking  changes  from  ^(j*) 
to  for  some  s*  e  S*  and  some  /,  /  »  1,2,..., AT.  Estimates  for  characteristics  of 

limiting  passage  times  can  be  based  on  measurement  of  passage  times  contained  in  these 
blocks.  Denote  the  state  space  of  the  process  {V(Tln):n  2  0}  by  Gl  and  set 

(4-8)  ~ 

/  -  1,2,..., AT. 

Denote  by  {Tn\n  2  1}  the  times  in  increasing  order. 

(4.9)  PROPOSITION.  Suppose  there  exists  S  >  0  such  that 

(4.10)  P{V(T\)  -  vjl  K(r‘_j) . nrj)}  2  8  a.s.. 


Then  P[V{Tn)  -  ^  i.o.}  -  1  for  all  /  -  1,2 . N. 


Arguments  analogous  to  those  given  in  Section  3  establish  Proposition  (4.9).  Using 
symmetry  of  the  spn,  the  idea  is  to  show  that 

(4.11)  P{V(T\)  =  vJmrJ.j)  »  •F(ri)  =  4<vo>} 


=  p{v&n)  =  J0\  n^)  =  ♦i(vi) . =  4<vo)} 


for  all  v0,v1,...,»n-1C  Gn.  (For  v  *  (s,s)  e  GN  we  write  <t»s(v )  =  (<f>5(i),<f>5(*))) 

Carry  out  the  simulation  of  a  0}  in  random  blocks  defined  by  the  successive 

random  times  i  0},  where 

(4.12)  fpk  -  inf  {7^  >  T'p  :  F(7^)  *  for  some  /,  /  *  1,2,.. .,JV}, 

k  2  1;  P0  -  0  and  7^  =  0.  Each  epoch  7^  corresponds  to  the  termination  of  a  passage 
time  with  no  other  passage  times  underway.  (Note  that  the  random  times  {7^:k  £  0}  do 
not  form  a  sequence  of  regeneration  points  for  the  process  {(2f(7)|),Pj|+l):/i  2  0}.) 


Set  ak  =  pk  —  Pk_v  k  i  1.  According  to  this  definition  ak  is  the  number  of  passage 
times  in  the  fcth  block.  Also  set 

“i 

rxv)  -  2  M'p 

j- 1 


and  denote  the  analogous  quantity  in  the  kth  block  by  Yk(f ),  k  z  1. 


(4.13)  PROPOSITION.  The  sequence  of  pairs  of  random  variables  {(Yk(f),ak):k  2  1}  are 
independent  and  identically  distributed. 

Proof:  As  in  the  proof  of  Proposition  (4.1),  observe  that  at  time  7^  defined  by  Equation 
(4.12)  a  passage  time  has  just  terminated  with  no  passage  times  underway  and  each  of  the 
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clocks  running  at  time  7^  +  was  set  or  can  be  viewed  as  having  been  probabilistically 
reset  at  time  7^.  Therefore  (Af(/):/  i  7^}  determines  the  distribution  of 
ajk+1  *  Pk+i  —  Pk  and  the  finite  dimensional  distributions  of  Ppk+i+ j  f°T  *  2  0-  The  joint 
distribution  of  the  clocks  set  or  reset  at  time  7^  depends  on  the  past  history  of  the  spn 
only  through  *(7^)  =  the  previous  marking  ^(j0),  and  the  trigger  transition 

It  follows  that  the  pairs  of  random  variables  {(.Yk(J),ak):k  z  1}  are  mutually 
independent. 

Recall  that  fn  is  the  time  of  the  nth  transition  firing  and  denote  by  e*  the  transition 
that  fires  at  time  n  2  0.  Also  recall  that  Cn  is  the  vector  of  clock  readings  at  time 
and  that  Cn  i  is  the  ith  coordinate  of  the  vector  Cn  for  E(Sn).  Let  Zj  z  0, 
Xj n€  S  and  CE  with  p(xk\xk_veit)  >0.  It  follows  from  the  definition  of  a 

symmetric  spn  that 

(4.14)  -  ♦><*„).  f.  s  v  #;  -  *«„_,)  - 

^n-l  S  Zn-V  en- 1  *  4<ej  } .  fl  *  4^1,)} 

1»~1  * 

-  'ii««  =  4<*,>-  t„  *  v -  4  <',.)•  *«„-,)  -  4<*„-.>. 

4  #„_i  “  4<',  ^ .  =■  4(»,,)} 

4—l  * 

for  all  /  ■  (Here  {- }  denotes  the  conditional  probability  associated  with 

starting  the  spn  with  marking  and  all  active  clocks  reset  at  time  /  =  0  according  to 

the  distributions 


PiC0j  Sc}*  7'(c;^(j0),eJ,^(j0),^£(e  )) 


for  c  i  0,  «,£  /*/{•}  denotes  the  corresponding  conditional  probability  when  the 

initial  marking  is  ^(jq).) 

Next  suppose  that  ,Y(0)  =  and  that  all  active  clocks  are  reset  at  time  t  -  0 
according  to  the  distributions  f(c;jQ,e|tJ0,e*),  e,e  E($q).  Set  AT1^)  =  ^(Af(/))  and 
x!(t)  «  $s(X (/)),  t  i  0.  Observe  that  for  each  sample  path  of  (AT(/):/  £  0}  and  all  niO, 

xl^«-i>€Ail  =  {*£>(*):  and  xhQeAp  =  {*£'(*):  s£A2] 

for  some  mi  if  and  only  if  Xl(Sn^{)  £  -d”7  an<*  X*(£n)  £  A™1  for  some  m,.  Similarly, 

x\!n_x)£B”'  =  {*£>(*):  jefij}  and  Af^C  -  {<$'(*):  *efi2} 

for  some  mj  if  and  only  if  €  B™1  and  X*(tn)  £  B2‘  for  some  m Since 

Sj(m)  =  inf  {f„  i  Tj(m ):  *(£„)£ 

and 

Tj(m)  -  inf  {f„  >  SV.jfm):  *(f„)  e  B%,  Jftf^)  e  S?} 

for  all  m.  Equation  (4.14)  implies  that 

(4.15)  />!{«!  -  n  +  1,  P’n+1  £  ^+1,  P'n  £  />;  s  yj 

-  ^{ctj  -  «  +  1,  /*;+,  *  yn+l ,  P'n  *  yn,...,  P\  s  yx} 

for  all  l  ■<  1,2,.. .,N,  +  l  2  and  «  *  0-  By  the  independence  argument  in  the 

first  part  of  the  proof  it  follows  that 
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(4.16)  *"  n  +  1.  ^fikkl  s  y*+\'  ^ Pkkl-\  S  *n . ^Pk+l  S  ?V  ~ 

-  f,{«,  - « + 1.  K+i s  s  >„ . /■',  s  >-,}  /■1{«iit)  -  4<>i)} 


for  all  n  2  0  and  /  =  1,2,...,#.  Using  Equation  (4. IS)  this  implies 


N 


2  ^*l{aJlc+l  *  w  +  ^t+1  *  yn+l'  *>pk+l-l  s  **Pk  + 1  s  -^l’  X^P^  “ 


■Y 


E  ^/{“l  -  «  +  1.  Cl  *  w  S  * . 'i  *  /I }  M*<^>  =  *W} 

;-i 


Y 


2  =  ft  +  1.  Cl  *  >„  +  1,  /»;  s  y^ . p\  *  yx}  P^Xitp  )  =  4(^)} 

/-I 


so  that 


^l{afc+l  *  n  +  1.  f3it+1  s  ->Wl*  ^^t  +  1_x  s  >n*"»  Ppk+l  S  -^l} 


/•1{a1  «  n  +  1,  Cl  s  yn+v  P'n  5  *„,...,  />i  s  yx] 


and  the  pairs  of  random  variables  {(T*(/),aA):*  z  1}  are  identically  distributed.  □ 


Standard  arguments  establish  a  ratio  formula  for  r{f)  =  £{/(/*)}. 


(4.17)  PROPOSITION.  Provided  that  £{tj}  <  «,  P[PeD{f)\  -  0  and  E{\f{P)  | }  <  », 


£l/TP)} 


E[YX(J)\ 
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With  these  results,  based  on  n  blocks  (cf.  Crane  and  Iglehart  [2])  a  strongly  consistent 
point  estimate  for  r{f)  is 


(4.18) 


r(n) 


Y(n ) 
a(n) 


and  an  asymptotic  100(1  -  2y)%  confidence  interval  is 


(4.19) 


*l_r  J(") 

5(n)  nl/2  ’ 


kn )  + 


S{n)  1 

a(n )  n1/2  J 


where  r2(n)  is  a  strongly  consistent  point  estimate  for  a2(f)  =  var  (Yx(f)  —  r(f)ax). 
Confidence  intervals  are  based  on  the  c.l.t. 


(4.20) 


nl/1  {r(/0-K/)} 
o<J)/E{ax\ 


WU) 


as  n-* oe. 


(4.21)  EXAMPLE.  In  the  token  ring  model  of  Example  (2.18),  consider  port  access 
times  measured  from  the  arrival  of  a  packet  for  transmission  by  some  port  until  the  start 
of  transmission  by  the  port.  This  sequence  of  passage  times  is  specified  by  the  four 
subsets  Ax  -  {(sv...jAN)e  S:  sx  =»  j3  -  0},  A2  «  {(jj,...,^)  €  S:  sx  -  1  and  *3  =  0}, 

B x  ■  {(j|,...^r^jy)  £  5":  Jj  *  0  and  r4jy  ™  1},  and  B2  3  {(■*!»... rf4jy)  £  S'.  Jj  “  1  and  =  0}. 
The  set  of  all  possible  markings  when  a  passage  time  terminates  or  is  not  underway  is 
H1  -  {(s1^2,...^4^)£S,:j4/_3s4/_1  -  1  or  r4/_2  *  1}.  Then  B\  #  0  and 
Jq  •  (1,0, 1,0,0, 1,0,0,.. .,0,1, 0,0) £  B2.  The  random  times  a  0}  correspond  to 

terminations  of  access  times  which  occur  when  there  is  no  packet  queued  for  transmission 
at  any  of  the  ports.  Propositions  (4.13)  and  (4.17)  hold  provided  that  the  packet 
interarrival  time  random  variables  are  exponentially  distributed.  (The  random  time  7^  is 


the  nth  time  at  which  port  /  +  1  observes  the  token,  n  a  0.  Note  that  {V(T Jn):n  a  0}  is  an 
irreducible,  finite  state  discrete  time  Markov  chain  so  that  PlViT^)  =  Vq  i.o.}  =  1.  It 
follows  that  P{V(Tn)  =  vj,  i.o.}  =  1  for  all  /  =  1,2 . N.) 

5.  STATISTICAL  EFFICIENCY 

Section  4  provides  two  estimation  procedures  for  passage  times  in  a  symmetric  spn. 
Each  of  these  procedures  rests  on  the  assumption  that  there  exist  e*  e  E,  sQeBv  and 
jq€/?2  satisfying  the  conditions  of  Proposition  (4.1).  The  regenerative  structure 
guarantees  that  Pn^P  as  n-»«>  and  the  goal  of  the  simulation  is  the  estimation  of 
r(f)  =  £{/!/’)},  where  /  is  a  real-valued  measurable  function.  (We  assume  that  the 
function  /  is  such  that  E{  \f\P)  | }  <  «  and  P{Pe  £>(/)}  =  0  so  that  ratio  formulas  for  iif) 
hold.) 

Estimates  for  r{f)  can  be  based  on  measurement  of  passage  times  {P\'.n  a  1}  and 
simulation  of  the  underlying  spn  in  regenerative  cycles  defined  by  the  times  rtn  at  which 
K(7^)  =  vq.  Alternatively,  exploiting  properties  of  a  symmetric  spn,  estimates  can  be 
based  on  measurement  of  passage  times  {P'„:n  2  1}  and  simulation  of  the  underlying  spn 
in  independent,  nonidentically  distributed  blocks  defined  by  the  times  at  which 
K(7^)e  {vq,...,v^}.  This  estimation  procedure  extracts  more  passage  time  information 
from  a  simulation  of  fixed  length  and  should  provide  estimates  for  r(/)  that  are  relatively 
more  accurate.  In  this  section  we  verify  that  this  is  indeed  the  case  by  showing  that  the 
resulting  confidence  intervals  are  shorter. 
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For  /  s  0  let  ml(/)  be  the  number  of  passage  times  {P\'-n  si}  completed  in  (0,/]  and 
denote  by  {/j£:/c  s  1}  the  indices  of  the  successive  termination  times  {T>rl:n  s  1}  at  which 
Vitn)  =  vj.  Set 

a\  =  -  m^T^i  ), 


r*(/>  =  £  A*}). 


y-m'd",  )  +  l 


Jk  i  1.  Also  set 


(a1^)  -  var  (yj{/)  -  K/)«J). 


Then  by  Lemma  (4.1)  of  Iglehart  and  Shedler  [5], 


V*/ 

'1/2(-r—  2  A*})-*/)) 

^(0  ' 

1  1/2  1/1 
(£{tJ})  aX(J)/E[a\ } 


M0,1) 


2  2 

as  f-*«  provided  that  £{(aj)  }  <  «  and  £{(t}(|/I))  }  <  *>.  Here  t*  =  T 'i  -  7^1  . 

Pk  Pk- 1 

Since  the  numerator  in  this  c.l.t.  and  the  limit  (^(0,1))  is  independent  of  the  transition 
)  and  the  markings  <^(j0)  an<*  $$(* (P  whi°h  define  the  cycles,  so  is  the  denominator; 
this  is  a  consequence  of  the  convergence  of  types  theorem  (Billingsley  [1],  Theorem  14.2). 


Thus,  the  quantity 


eHf)  =  (£{rj})  VG0/£{aJ} 


is  an  appropriate  measure  of  the  statistical  efficiency  of  the  estimation  procedure  based 


on  cycles. 
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blocks. 
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Note  that  when  the  passage  times  {P\:n  a  1}  are  used  to  construct  point  and 
interval  estimates  for  r(J),  the  half-length  of  the  confidence  interval  is  proportional  to 
ex(f),  and  when  the  passage  times  {P^rn  ^  1}  are  used,  (with  the  same  constant  of 
proportionality)  the  half-length  of  the  confidence  interval  is  proportional  to  e(f). 
Proposition  (S.3)  asserts  that  under  mild  regularity  conditions  on  the  function  /, 
eif)  *  eX(J). 

(5.3)  PROPOSITION.  For  all  functions  /  such  that  E{  \J\P)  1}  <  »  and  P[Pe  D(J)\  =  0, 
eif)  £  ex(J). 

Proof:  It  is  sufficient  to  show  that 

(5.4)  (<r(/))2  S  N2  (a^/))2 
and 

(5.5)  E{ai]  -  NE{a\}. 

To  establish  Equation  (5.4),  for  t  2  0  set 

m(/) 

m0  =  2  /IP)  -  rif)m(t). 

j- 1 

Now  observe  that  {*(/):/  £  0}  is  a  regenerative  process  by  Proposition  (2.19)  and  that, 
with  respect  to  this  process,  [W(0’t  a  0}  is  a  cumulative  process  in  the  sense  of 
Smith  [17]  with 

E{W(fpl)  -  Wfpi  )}  =  E[Yk<J)  -  r{f)ak\  =  0. 


Thus,  by  Theorem  S  of  [17], 


(5.6)  am  ,-ar-^(,))  -  «a£. 

'  Elrjl 

Next  recall  that  {P'n’-n  a  1}  is  the  sequence  of  passage  times  P\,P\,...,P^ ,p\ 
enumerated  in  termination  order  and  therefore 


«(f)  N  f m\t ) 

2  KP'P  -  «  2  <  2  KP*)  -  rif)m\C) 

j-1  U\  lj-l 


where  ml(t )  is  the  number  of  passage  times  [P^in  a  1}  completed  in  (0,/].  Now  set 


«'(/) 

rf(i)  -  2  /*/>£  -  r{f)m\i) 
j- 1 


so  that 


N 

mt)  =  5;  ^(/) 

/-i 


and  by  the  Cauchy-Schwarz  inequality 


JN  1/2 

var  (1^(0)  *  2  var  (W^(/))  +  2)  {var  ( W*Q))  var  (»**(/))} 

/-l  j+l 


S  1/2 

2  {var  (W^(/))} 

/- 1 


lim 

f-*«0 


var  (W^(/)) 

t 


E{  tJ] 


Equation  (5.4)  follows  since 


for  /  =  1,2,...,M  To  see  this,  fix  /  and  let  *  1}  be  the  indices  of  the  successive 

termination  times  [7^:n  *  1}  at  which  V(?n)  *  vq.  Observe  that  is  a 

cumulative  process  so  that  by  Theorem  8  of  [17], 


lim 

r—*o 


var  (W^(/)) 
/ 


(g'Cfl)2 
£{t^}  * 


where  rk 


'ffj  —  T'oi  and 
p*  pk-  i 


(a^/))2  =  var  (l^CO  -  rC04), 


with  alk  »  ml(T,fj)  —  nt'ct-i  )  and 

Pk  Pk“! 

>*</>  -  s  /vj). 

)  +  l 

*“i 


The  definition  of  a  symmetric  spn  implies  that  E{t12\  -  E{t\]  and 
(*'(/))2  =  var 


To  establish  Equation  (S.S)  set 


N 

m(t)  »  £  m\t) 
/-I 


and  observe  that  [m7(/):/  iOj  and  {/*(/):/  2  0]  are  cumulative  processes  with  respect  to 


{Af(/):r  a  0}.  Moreover, 


Again,  since  the  spn  is  symmetric,  E{ Tj}  —  E{t\)  and  £{aj}  =  E[a\\  so  that 


lim 

f-»«0 


t 


£{«}} 

E{r\] ’ 


/  -  1,2 □ 
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Figure  1.  Token  ring. 


Figure  2.  SPN  representation  of  token  ring. 
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